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(~| I CFHTLS optical photometry has been used to study the galaxy luminosity functions of 

I' 14 X-ray selected clusters from the XMM-LSS survey. These are mostly groups and poor 

clusters, with masses (M500) in the range 0.6 to 19 x IO^^Mq and redshifts 0.05 < z < 0.61. 
Hence these are some of the highest redshift X-ray selected groups to have been studied. 
Lower and upper colour cuts were used to determine cluster members. We derive individual 
^ I luminosity functions (LFs) for all clusters as well as redshift-stacked and temperature-stacked 

LFs in three filters, g', r' and z', down to M — —14.5. All LFs were fitted by Schechter 
functions which constrained the faint-end slope, a, but did not always fit well to the bright 
end. Derived values of a ranged from —1.03 to as steep as —2.1. We find no evidence for 
upturns at faint magnitudes. Evolution in a was apparent in all bands: it becomes shallower 
with increasing redshift; for example, in the z' band it flattened from -1.75 at low redshift to 
QQ ', -1.22 in the redshift range z =0.43-0.61. Eight of our systems lie at z ^ 0.3, and we combine 

■ these to generate a galaxy LF in three colours for X-ray selected groups and poor clusters 

at redshift 0.3. We find that at z ^ 0.3 a is steeper (-1.67) in the green (g') band than it is 
(-1.30) in the red (z') band. This colour trend disappears at low redshift, which we attribute 
to reddening of faint blue galaxies from z ~ 0.3 to z ~ 0. We also calculated the total optical 
I luminosity and found it to correlate strongly with X-ray luminosity (Lx oc Lpprp), and also 

with ICM temperature {Lqpt « T^^^^), consistent with expectations for self-similar clusters 
with constant mass-to-light ratio. We did not find any convincing correlation of Schechter 
parameters with mean cluster temperature. 
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1 INTRODUCTION 

Most of our knowledge of galaxies is based on observations of 
the local universe, although distant universe observations have also 
provided a wealth of information. Statistical studies of galaxies 
at high redshift are mostly limited to rich galaxy clusters mainly 
due to observational limitations. Galaxy clusters are important cos- 
mological environments where key galaxy transformation such as 
stripping and strangulation occur. However, in the hierarchical for- 
mation of structure rich clusters are the latest structures to be 
formed. Lower mass systems or galaxy groups may have been the 
place where galaxies experience a substantial degree of evolution 
through processes such as mergers and tidal interaction, as a re- 
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suit of the higher efficiency of these processes in the lower velocity 
dispersion environment of groups. 

The Galaxy luminosity function (LF) - the number of galaxies 
per unit volume in the luminosity interval L to L + dL ~ has been 
widely used to study the formation of galaxies and the evolution 
of galaxy populations with redshift. It is also an excellent statis- 
tical tool for describing how different environments influence the 
properties of galaxies. 



Both the bri ght end (iBower et alj ( l2006h, iNaab et alj ilOOH) ) 
and the faint end jMarzke et alj ( ll994h . TKhochfar et alj ilOOHi ) of 
the LF have been the subject of in-depth studies, as they offer strong 
observational constraints for models of galaxy formation and evo- 
lution . While the bright end of the LF is affected by AGN feed- 
back teower et alj i2006l) ), the faint-end slope is predominan tly in- 
fluenced by feedback from supernovae toekeletalj ( ll986l) ). and 
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provide a direct indicator of thie significance of dwarf galaxies, 
which are expected to behave differently in rich and poor clusters. 
Multi-colour LFs, in particular, probe the history of the faint galaxy 
population, includi ng its star formation history - see for example, 
lAdamietalJ ( l2007l) . 

The vast majority of studies of the galaxy LF give faint-end 
slopes in the range ~ — 1 to ~ —2. Most of these have lim- 
ited magnitude depth [M > —16) and recent de ep studies are 
mostly confined to rich l ocal clusters (See Table 1 in lPopesso et alj 
( l2005ij) and Table A. 1 in lSoue et al] j2008l) and references therein). 
These studies not only disagree on the value of the faint-end slope, 
but the y also disagree on th e exact form of it, as some stud- 
ies (e.g. iGonzMez et al] ( |2006|) ) found upturns; a single Schechter 
function was not an adequate fit to the faint end, and a double 
Schechter function was required to give a reasonable fit. The ex- 
istence of these upturns is very sensitive to the method used to de- 
termine galaxy membership, with some approaches including spu- 
rious galaxies or excluding genuine cluster members due to their 
low surface brightness. 

The evolution of the faint-end slope is hard to study, mainly 
because the number of faint galax ies detected decreases sharply 
with increasing redshift. iLiu et al. 

I I2OO8) found that the faint-end 
slope of a field galaxy population became shallower with increas- 
ing redshift (up to 2; = 0.5) for all galaxy spectral types. However, 
to account for the photometric redshift errors of the galaxies, they 
weighted the galaxies as probability-smoothed luminosity distribu- 
tion at the redshift at which they were measured. This places an im- 
portant caveat on the interpretation of their data, and hence on their 
results. On the other hand, simulations by iKhochfar et al. 
show a measurable dependence of the faint-end slope of the galaxy 
luminosity function on redshift. However, most of this dependence 
is seen over a relatively large redshift range, Az > 2. Furthermore, 
it is hard to discriminate galaxy environments in such studies. 

X-ray surveys remain one of the most popular methods of 
finding galaxy systems. Due to the strong density dependence of 
X-ray emissivity. X-ray cluster selection is much less vulnera- 
ble to contamination along the line-of-sight t han optical method s. 
The XMM-Large Scale Survey (XMM-LSS) jPierre et alj ( l2004h ). 
a contiguous X-ray survey, has a well-defined selection function 
which is used to produce a sample of galaxy groups to study their 
intra cluster medium and galaxy properties at medium to high red- 
shift. IPacauTeTd] ( |2007|) have presented a study of a sample of 29 
galaxy systems from the XMM-LSS survey, drawn from an area of 
5 deg^ out to a redshift of z = 1.05. The cluster distribution peaks 
around z = 0.3 and T=1.5 keV, half of the objects being groups with 
a temperature below 2 keV. 

In this paper, we use the XMM-LSS optical follow-up obser- 
vations to study the evolution of the galaxy luminosity function 
in galaxy groups and poor clusters since z ~ 0.6. Given the ob- 
servational biases - distant groups are more massive and hotter - 
we study whether the redshift dependencies are weaker or stronger 
when the intrinsic properties of the systems, for instance, intraclus- 
ter medium temperature, are taken into account. 

By using a deep (m^/ = 24) optical survey of X-ray selected 
galaxy clusters up to redshift of z = 0.61, we aim in this paper 
to clarify the debate on the faint-end slope of the LFs of low-mass 
(M500 < 20 X 10^'^ A/0) galaxy clusters (or groups), and to explore 
the existence of any dips, or upturns at the faint end, and to estab- 
lish whether the slope shows trends with redshift or intracluster 
medium temperature. Comparison with previous results can help to 
elucidate the universality of galaxy cluster LFs. Furthermore, the 
scaling relation of total optical luminosity with temperature and 



X-ray luminosity for our cluster sample can shed light on the mass- 
to-light ratios of low-mass systems when compared to rich clusters. 

The paper is constructed as follows: In section 2, we describe 
the optical catalogue used to calculate the LFs. Then, we describe 
the data reduction and the method used to construct the colour- 
magnitude diagrams (CMD) and the subsequent LFs, and the tech- 
nique adopted for the background subtraction. In section 3, we de- 
scribe our results, starting with the individual cluster LFs, and then 
the redshift- stacked clusters and temperature- stacked clusters. In 
section 4, we discuss our results and compare them with other stud- 
ies. Finally, in section 5, we summarise our conclusions. 

Throughout this a rticle, we ad o pt the cosmological pa- 
rameters estimated by ISpergel et alj hoOTi) . namely: Ho = 
73 km s-i Mpc-\ = 0.24, Ha = 0.76. 



2 DATA 

2.1 Observations 

Optical photometry of the XMM-LSS survey was obtained from 
the Canada-France-Hawaii Telescope Wide Synoptic Legacy Sur- 
vejQ, referred to as the CFHTLS Wide survey. Data were obtained 
in five passbands (u* , g' , r' ,i' , z') down to a nominal magnitude 
limit of i' = 24.5. Of the 19 deg^ of CFHTLS Wide data avail- 
able in the WI survey area, 4 deg^ overlap with th e X-ray selected 
cluster catalogue presented bv lPacaud et aljf2007l) . Hence our pho- 
tometric data are drawn from four 1° x 1° catalogues derived from 
the survey data. 

The data used in this paper are based upon the reduction proce- 
dure outlined in Hoekstra et al. ( 20 06). Source extraction and pho- 
tometry were performed using SExtractor v2 . 5 . (Bertin 
and Arnouts 1996). Zero point information for sources detected 
in the CFHTLS Wide field survey WI area was extrapolated from 
common sources detected in the Sloan Digital Sky Survey equato- 
rial patch which overlaps the southern edge of the WI area. 

XMM-LSS Class 1 (CI) clusters are a well-controlled X-ray 
selected and spectroscopically confirmed cluster sample. The cri- 
teria used to construct the sample guarantee negligible contamina- 
tion by point-like sources. The observations of the clusters were 
performed in a homogeneou s way (10-20 ks exp osures). For full 
details of the CI sample, see lPacaud etalj j2007l) . The mam prop- 
erties of the sample are shown in Table [T] Detailed information on 
the CI selection process can be found in Pacaud et al. (2006). 17 
out of the 29 XMM-LSS CI clusters are covered by the CFHTLS 
Wide field survey. In this paper, we study the luminosity functions 
of 14 of these 17 clusters - dropping the three with the highest 
redshifts (clusters with XLSSC numbers 2,29, and I) because their 
photometric data is too poor to allow useful constraints to be ob- 
tained. 

2.2 Analysis 

Galaxies were detected by SExtractor dBertin & AmoutsI h996t) ). 
Luminosity functions (LFs) were produced in three of the five 
CFHTLS (u* , g', r' z') filter bands, namely, g', r' and z' . To 
determine the completeness of the LFs, we took into account the 
limiting apparent magnitude in each field. The completeness limits 
for each filter was determined using the apparent magnitude LFs 
of all data (down to the faintest magnitudes available) for each C 1 

1 See |http://www.cfht.hawaii.edu/Science/CFHLS7] 
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cluster individually. Variations in seeing and exposure time across 
the CFHTLS fields used here are small, and it was found that for 
each filter there was a common completeness limit at which all LFs 
started to drop below the faint end power law slope. Note that the 
LF turn-up reported by some authors (see section 4.2), which could 
potentially introduce an error into this method for estimating com- 
pleteness, falls beyond the faint limit of our LFs (e.g. at -16 in g' 
band), except in our three closest clusters, and hence cannot seri- 
ously bias our estimates of the completeness limits. The complete- 
ness threshold magnitudes for the three filters g', r' and z' were 
found to be 24, 23.5 and 23, respectively. These values are also 
consistent with results based on comparison of the number counts 
per field to deeper data from the CFHT LS Deep Field and CCCP 
Megacam observations dUrquhart et al] |2009)). 

Each entry in the catalogues is associated with a FLAG value 
which indicates the degree of reliability of the data. Flag is a short 
integer, and a value of denotes good data. The more unreliable the 
data is, the higher the FLAG value becomes. We included all cat- 
alogue entries with FLAG < 3, which includes sources with very 
close and bright neighbours or some bad pixels and sources which 
are originally blended with other sources. This may admit some 
problematic galaxies but this is better than excluding many genuine 
cluster members, because many clusters contain significant number 
of blended sources. Factors that may raise the FLAG to > 3 include 
sources with saturated pixels, truncated sources, incomplete or cor- 
rupted data and sources with memory overflow during deblending 
or extraction. Catalogue entries with FLAG > 3 constitute only 
~ 5% of the total number of entries, and were all removed. Many 
of the removed entries are fainter than the threshold magnitude and 
hence would have been removed anyway. 

Each entry in the catalogue also includes a stellarity class 
value, STAR, with values ranging from to 1 . The lower its value, 
the more likely the detected object is a galaxy. Data points with dif- 
ferent STAR values were checked by IRAF and their radial profiles 
were examined to see if they matched the typical profile of a star 
or a galaxy. Typically, objects with STAR > 0.85 were found to be 
stars, whilst those with < 0.85 were galaxies. Therefore, only cata- 
logue entries with STAR values of 0.85 or less were included when 
constructing the LFs. 

Spectral temperatures of the XMM-LSS clusters, and the re- 
sulting values of -R500 (which allow for the evolution i n critical den- 
sity w ith redshift), were measured by the Saclay team jPacaud et al.l 
( 120071) ). To construct colour-magnitude diagrams, we selected all 
galaxies within a circle of radius, R* = 1.5 x -R500 of the clus- 
ters. This radius limit, R* , represents an estimate of i?200- Colour- 
magnitude diagrams were produced for all 14 clusters for colour 
bands: g', r' and z' . The factor, 1.5 does not have a large effect on 
the fitted parameters of the Schechter function; we compared the 
results of 1.0 x -R500 to 1.8 x -R500 and found that the faint-end 
slope, a, was only changed within its 1 a errors. 

The CMD were used to colour select galaxies which might be 
cluster members, hence reducing the background due to interlopers. 
The colours used for this were U2 — g2 versus g^ron, g'2 — r'2 versus 
rJ,jQn, and i'2 — z'2 versus z'^^^^-^ for the three filters respectively, 
where the subscript 2 refers to the 2 arcsec aperture used in the 
magnitude measurements. To define and select cluster members in 
the CMD, we defined upper and lower colour cuts and only galaxies 
between these two lines were used to produce the LF, as galaxies 
outside these limits were most likely not cluster members. To define 
these two colour cuts, we first defined the red sequence line in the 
CMD and then pushed this line up and down to allow for statistical 
errors, and for the likely range of galaxy colours. 



To define the red sequence line, we first defined the slope and 
then its Y-intercept. We checked that the BCGs lay at the centre 
of the X-ray emission in all clusters, and then calculated the red 
sequence slope in each case by fitting a straight line to the bright 
galaxies in the CMD. Bright galaxies are defined as those with mag- 
nitude ranging from that of the brightest cluster galaxy, ttibcg, to a 
magnitude of mBCG + 3, inclusive. We found that the slope of the 
red sequence line for the 14 CI clusters showed a mild trend with 
redshift: it was s teeper for high - redsh ift clusters. A similar trend 
was observed by iGilbank et alj ( l2008h and attributed to a deficit 
of faint red galaxies at high redshifts, consistent with the galactic 
downsizing picture. 

We divided our cluster sample into two redshift ranges: low- 
redshift clusters (z < 0.2) and intermediate-redshift clusters 
(0.2 < z < 0.61). Clusters from each group share a common red 
sequence line slope with a small variation. The common slope for 
the low-redshift range was -0.007 and -0.025 for the second range. 
Instead of using a different red sequence slope for each cluster, we 
used the common slope of the redshift groups for all clusters be- 
longing to that redshift group. 

The Y-intercepts of the red sequence lines were different for 
each galaxy cluster and depended on the average colour of the 
bright galaxies as defined above. To fix the value of the inter- 
cept for each cluster, the red sequence line was normalised so as 
to pass through the point in the CMD which has a magnitude of 
m-BCG + 1-5 and colour equal to the average colours of the bright 
galaxies. This point and the value of the slope completes the defi- 
nition of the red sequence line. 

Both upper and lower colour cuts have the same slope as the 
red sequence line. In order to define the upper colour cut, we have 
to determine the upper (red) limit to the cluster red sequence. We 
took into account the statistical scatter of the colours of the faintest 
galaxies on the red sequence. These galaxies are defined as those 
inside a 1.0 x 0.1 (magnitude by colour units) box in the CMD 
centred on the faint end of the red sequence line (see figure[T]). The 
size of this box was chosen to include the faintest galaxies most 
probably belonging to the red sequence after studying the CMD of 
the CI sample. The expected scatter of these galaxies, a, is calcu- 
lated by averaging their colour errors, that is the Y-axis errors in the 
CMD. The upper colour cut, is then taken to be the red sequence 
line pushed upward by 2a. By taking into account this scatter, we 
ensure that almost all genuine cluster red sequence galaxies should 
fall beneath the red cut, since the statistical error on the brighter 
galaxies will be smaller. 

Similarly, the lower (blue) colour cut is the red sequence line 
pushed downward in the CMD. In this case, the shift has to ac- 
count for both statistical scatter, and for the fact that late-type clus- 
ter galaxies are intrinsically bluer than red sequence galaxies. The 
shift was therefore taken to be —{2a + A). Where A is the the- 
oretical colour difference between ellipticals and spirals. This was 
estimated using a simple model which calculates what colour late- 
type galaxies would have when redshifted by different amounts, 
as described in King & Ellis (1985). A is a function of redshift 
only and the redshift of the galaxy cluster was used to determine 
its value. This method of estimating A ignores any intrinsic evo- 
lution in the colour offset between red and blue cluster gala xies. 
However, the detailed COMBO- 17 study of iBell et alj j2004 l) (see 
their Figure I) shows that the colour difference between blue and 
red sequence cluster galaxies changes little over the redshift range 
(0-0.7) spanned by our clusters. Figure[T] shows an example of our 
use of colour cuts for selection of cluster galaxies. 

Of course, background and foreground galaxies will still con- 



4 A. Alshino et al. 




Figure 1. Colour-magnitude diagram of cluster 25 (redshift=0.26). All 
galaxies (crosses) on the left of the vertical dotted Une are the bright galax- 
ies with magnitude < rriBCG + 3. The red sequence line (the solid line) 
is defined by the point (diamond) with magnitude of mecG + 1-5 and 
colour equal to the average colours of all bright galaxies and by the slope 
of -0.025. The statistical scatter, cr, is estimated by the average colour er- 
rors of the galaxies within the 1.0 X 0.1 box on the faint end of the red 
sequence line. In the case of cluster 25, cr = 0.07. The dashed lines are 
the upper and lower colour cuts. The upper colour cut is the red sequence 
line pushed upward by 2a and the lower colour cut is the red sequence line 
pushed downward by 2a + A, where A = 0.175 for redshift of 0.26. See 
text for definition of A. 



taminate the sample after the colour cut has been applied, and this 
contamination must be estimated and removed statistically. For this 
purpose we used all data in the catalogue to which a given galaxy 
cluster belonged. In addition to simple Poisson fluctuations, un- 
certainties in removing background and foreground galaxies arise 
from large scale structure. To quantify the extra fluctuations arising 
from this, we proceeded as follows. The whole catalogue 1° x 1° 
area was divided into smaller blocks with areas comparable to that 
of the cluster in question. Any of these blocks covered mostly (60% 
or more, by area) by a galaxy cluster, were considered to be dom- 
inated by a cluster and hence were discarded from the background 
calculation. Blocks covered by clusters to an extent less than 60%, 
were not discarded but the portion covered by the R* circle of any 
galaxy cluster was removed, so the final blocks used have some- 
what different areas. 

For each background block, an LF was produced in just the 
same way as for the cluster itself. The same values of the upper and 
lower red sequence limits of the galaxy cluster in question, were 
applied to all its background block areas, so galaxies beyond those 
limits were removed. The application of colour cuts to both source 
and background fields reduces the noise level in both of them, and 
hence in the final background-subtracted LF. 

We then divided each background block LF by its area, added 
them and normalised the resulting single LF to the area of the 
galaxy cluster in question. The error bars on the averaged back- 
ground LF were calculated from the scatter of the individual block 
LFs contributing to it. This method of estimating the background 
has some advantages over the more conventional background esti- 
mation method using an outer annulus around the galaxy cluster, 
since it uses a large background region, and the error estimate al- 
lows for the variance arising from the large scale structure. Finally, 
for each cluster we subtracted its composite background LF from 
the cluster LF, and propagated the errors. 



Apparent magnitudes were converted to absolute magnitudes, 
using the distance for each cluster, and applyi ng K-corrections cal - 
culated from Table 3 (for Hubble type E) in lFrei & Giiiuil ( Il994l) . 
The use of early-type K-corrections is common in cluster studies, 
and justified by the dominance of early-type galaxies in clusters. 
However, if there were a systematic trend in early-type fraction 
with magnitude, then this could lead to some distortion of the LF 
slope. To quanti fy the maximum po ssible effect, we note that, us- 
ing the tables in lFrei & GunnI ( 1 19941) . the K-correction for Hubble 
type E at z=0.6 for z' is 0.37 while it is 0.05 for Hubble type Im). 
Assuming (very conservatively) a systematic change from a 100% 
early-type to 100% late-type population across the faint end slope 
of our LFs, the impact of a differential error of 0.3 magnitudes on 
our determination of a would still only amount to Aa ~ 0.04, 
which is small compared to the trends in alpha wh i ch rep resent 
some of our main results. The tables in lFrei & GunnI ( 1 19941) apply 
to SDSS filters, which differ slightly from the corresponding Mega- 
Cam filters. The resulting differences in K-corrections are much 
smaller than the differences between early-type and late-type galax- 
ies (about 0.03 at z=0.1 and 0.06 at z=0.6), and will have negligible 
effect on our derived LF slopes. 

Finally, the data were binned into bins of width 0.5 magnitude 
(experiments showed that this bin size was a good choice in terms 
of fit quality and parameter confidence regions), and the result- 
ing LF s were then fitted by a Schechter function model jSchechte3 
Il97i)), 

(j>{M)dM = 0.4/n(10)(^*e"^X^+"dM, 

where X = iq-O-Hm-m") ^ characteristic magnitude, 

(j)* is the c harac teristic number density and a is the faint-end slope, 
lLinetalJ lll996). Contour plots of the Icr, 2a and 3cr confidence 
levels of a and M* were also produced. The errors in the text and 
tables refer to the la errors. We also calculated the total optical 
luminosity Lqpt of each cluster by integrating the fitted Schechter 
function from 5 x M* to -16. 

In addition to single LFs for each galaxy cluster in our sam- 
ple, we produced stacked luminosity functions. The radius used to 
determine the volume is the R* of the cluster. Before stacking dif- 
ferent clusters together, to correct for the evolution in the critical 
density of the universe, we multiply the LF by 

Pc[z = 0) 

Pc{z = Zcl) ' 

where z is the redshift, pc is the critical density of the universe, a 
function of z, and Zd is the redshift of the cluster. This correction 
is necessary for high-redshift clusters if stacked with low-redshift 
clusters to scale the galaxy density in each cluster to the density at 
redshift=0. The faintest magnitude bin is not necessarily the same 
for each cluster and to account for this, we divided the total num- 
ber of galaxies in each magnitude bin by the summed volume of 
galaxy clusters that contributed to that bin only. The stacked LFs 
should enable us to study the evolution of the LF with redshift and 
to explore any differences between clusters of different tempera- 
ture. 



3 RESULTS 

3.1 Individual cluster luminosity functions 

The fitted values of the three Schechter parameters a and M* and 
Lo PT of the individual C 1 clusters are presented in Table |2| The 
LF plots with the associated la, 2a and 3(7 contours in the A/* -a 
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Table 1. List of the 17 CI galaxy clusters covered by CFHTLS optical survey and their properties sorted according to their redshifts jPacaud et alj i2007h ). 
The three highest redshift clusters (2,29 and 2) though covered by the survey, were not included in our analysis because their data were too poor to yield useful 
fits. 



plane for passbands r' and z' are shown in Figures |AT] and |BT] re- 
spectively. For some of the CI clusters, the fitting failed to constrain 
some of the parameters, M* in particular, due to poor statistics or 
the lack of any well-defined turnover in the LF at the bright end. 
For these clusters the LF and best fit are presented without any ac- 
companying confidence contour plot. These LFs are placed at the 
bottom of the figures. For clusters with unconstrained M* , Lqpt 
was also not constrained, because its value depends on both a and 
M* . Therefore, we excluded these clusters in the part of the analy- 
sis related to Lo pt- 

The average values of ct for our sample of clusters are 

-1.70 ± 0.10, -1.64 ± 0.04 and - 1.43 ± 0.03 

for the g' , r' and z' passbands respectively. The correlations be- 
tween Loft, ol and M* and re dshift, temperature (T ) and the X- 
ray luminosity {Lx) taken from lPacaud et al.| teOOT^), were tested 
using Pearson's correlation coefficient. These coefficients are com- 
puted from the ratio of the covariance of the tested variables, X 
and Y, to the square root of the product of the variances of these 
variables, i.e. 

_ GOV {XX) 

■sJVAR{X) * VAR{Y) ' 

This correlation coefficient measures the linear correlation, if 
it is 1 or -1 then the two variables are perfectly positively or neg- 
atively linearly correlated, respectively. To compute the upper and 
lower 1 G errors on the correlation coefficient r, we used Fisher's Z 
transformation: Z = tan^^r. The strongest correlations found are 
those between Lqpt and T and between Lqpt and Lx, both of 
which are expected from the scaling relations of galaxy clusters. In 
our sample, they both have a correlation coefficient of at least 0.9, 
see Table |3] 

Because higher redshift clusters are more difficult to detect 
than nearby ones, they will tend to be more m assive and hence hot - 
ter than typical nearby clusters, see Figure 3 in lPacaud et ai] l l2007h . 
This (Malmquist) selection effect is present in any deep cluster sur- 
vey. To account for the T ~ z correlation arising from this selection 



effect in the CI sample, for each correlation coefficient between a 
quantity and T or z, we have also calculated the partial correla- 
tion coefficient between the same two quantities, which attempts to 
remove any part of the correlation which arises due to the intrinsic 
trend in T with z within our sample. For this we used an Interactive 
Data Language (IDL) routine, p. correlate. pro to compute the par- 
tial correlation coefficient. This uses the following method, which 
to be concrete we explain using the example of the correlation be- 
tween a and redshift. Let a and redshift z are the variables of pri- 
mary interest, whilst temperature T is a third variable whose effects 
we wish to remove. First, the routine calculates the residuals after 
regressing a on T; these are the parts of a that cannot be predicted 
by T. Likewise, it calculates the residuals after regressing z on T. 
Finally, the partial correlation coefficient between a and z, adjusted 
for T, is the correlation between these two sets of residuals. 

The results of our correlation analysis for the unstacked clus- 
ters, are tabulated in the top section of Table[3l The correlation co- 
efficients between the faint-end slope, a of the individual clusters 
and redshift are 0.44±0.27 for the r' band and 0.54±0.25 for the z' 
band. These coefficients, including the coefficient for the g' band, 
get stronger after the application of the partial analysis and the er- 
rors on the coefficients become smaller. This strongly suggests evo- 
lution of Q with redshift in our sample. We will further scrutinise 
this possibility in the section of redshift-stacked clusters, because 
stacking LFs of clusters with similar redshifts should lower scatter 
in the data and provide a means to probe possible trends. M* also 
shows a negative correlation with redshift and with temperature but 
these correlations become insignificant in a partial correlation anal- 
ysis. 



3.2 Global scaling relations 

The relationships between the global cluster properties, Lo pt, Lx 
and T provide a probe of cluster self-similarity. Lqpt is strongly 
correlated to the temperature of our clusters - the correlation co- 
efficients between Lqpt and T are 0.95±0.06, 0.96±0.04 and 
0.97±0.03 for the g', r' and z' bands respectively, whilst the partial 
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Figure 2. Correlation diagrams of Lqpt versus X-ray gas temperature, 
T (top panel) and Lqpt versus X-ray luminosity, Lx (bottom panel) of 
CI clusters for passbands g' (stars), r' (triangles) and z' (squares). Clusters 
with unconstrained M* and hence unconstrained Lqpt were excluded. 



correlation coefficients for the same quantities, factoring out the ef- 
fects of z, are 0.87±0.16, 0.89±0.11 and 0.92±0.06, see third row 
in Table [3l The removal of the z effects has lowered the values of 
the coefficients but they are still high and significant. Correlation 
between Lo pt and Lx is also quite strong: 0.92±0. 1 1 (g' band), 
0.93±0.07 (r band) and 0.90±0.08 (z' band). 

In Figure |2l we plot Lqpt versus T (top panel) and Lqpt 
versus Lx (bottom panel). We calculate the slopes for the se 
plots using the Fortran package ODRPACK (' Akritas et al.| ( 1 19961) ). 
which uses numerical orthogonal distance regression method to 
minimise perpendicular distances between points and the fitted line. 
One advantage of this is that the slope value will not change if the 
quantities in question switch axes. In addition, ODRPACK takes 
into account errors on both X- values and Y- values which are avail- 
able for Lo PT, T and Lx ■ 

The logarithmic slopes for the Lqpt — T relation for the 
three filters g' , r' and z' , respectively are 1.57 ± 0.17, 1.51 ± 
0.17 and 1.79 ± 0.12, giving an average value of 1.62 ± 0.11. 
For the Lqpt - Lx relation, the slopes are 0.47 ± 0.07, 0.43 ± 



0.08 and 0.50 ± 0.07, and the average value is 0.47 ± 0.05. Note 
that the slopes do not differ significantly for the three filters, ex- 
cept the slope of Lo pt versus T in the z' filter. Such relations be- 
tween Lq pt on one hand, and Lx and the gas temperature on the 
other, are expected because richer and hence more luminous clus- 
ters have deeper gravitational potential wells which in turn raise 
the ICM temperature and its X-ray output by adiabatic compres- 
sion and shocks generated by supersonic motion. We will discuss 
this further in Sec. 4.7. 

The correlation coefficients between Lqpt and redshift are 
high (all above 0.8), but when the effects of temperature are re- 
moved they become insignificant in at least two of the filter set, 
therefore, this correlation is most likely due to selection effects, 
Malmquist ejfect, and does not reflect any genuine relationship be- 
tween LoPT and z. Correlations between a and M* and T, z and 
Lx were also computed, but none of those showed significantly 
high values. 

Following the above analysis of trends in the properties of in- 
dividual clusters, we now perform a stacking analysis, grouping 
clusters first by redshift, and then by temperature. This provides 
LFs of higher statistical quality, enabling the behaviour to be ex- 
amined in greater detail. 

3.3 Redshift-stacked clusters 

The 14 CI clusters span a redshift range 0.05 to 0.61. This range 
was divided into five redshift bins: 0.05-0.14, 0.26-0.26, 0.29-0.29, 
0.30-0.32 and 0.43-0.61. The number of clusters in each bin ranges 
from two to four. The redshift ranges of these bins were chosen 
according to two criteria: first, the redshift range of the combined 
clusters was not too large, and second we required adequate data 
quality in each bin, to allow a well-constrained Schechter function 
fit. We kept the number of bins to at least five because a smaller 
number of bins increases the errors on the correlation coefficients. 
Plots of the redshift-stacked data with fitted Schechter functions 
for the three photometric bands are shown in Figures [TO] [TT] and 
II2I and results of the fits are given in Table|4] 

The faint-end slope, a of the Schechter function of the stacked 
data shows an evolutionary trend, becoming less steep with increas- 
ing redshift. Three of the redshift bins (0.26-0.26, 0.29-0.29, 0.30- 
0.32) have very similar redshifts and in general the a values for 
these three bins agree within their errors. 

The Pearson and partial correlation coefficients were calcu- 
lated for a and z, see Table [3] The coefficients are high (> 0.88) 
but with relatively large errors, mainly due to the small number of 
bins. The partial correlation analysis lowered the values of the coef- 
ficients and enlarged the errors. Evidence for evolution in a is seen 
in all three bands, arising primarily from the fact that the faint-end 
slope is steeper (a — —1.75 to -1.8) in the low z bin than in the 
higher redshift bins. 

One obvious concern in probing evolutionary trends in the 
Schechter function fits is that the fitted magnitude range decreases 
systematically with redshift, due to the apparent magnitude limit of 
our data. A second effect which might bias a is that within a given 
redshift bin, the contributing clusters are probed down to different 
absolute magnitudes, according to their distance. Hence at the faint 
end, clusters may progressively drop out of the stacked LF. This is 
especially the case for the lowest and highest redshift bin, which 
are both much broader that the three bins at z ~ 0.3. To show the 
scale of this latter effect, we have drawn a vertical dotted line on 
each of the stacked LF plots l|9][lO][II][l2l[T3l[l4]and[T5) to show 
the faintest magnitude to which all clusters in the bin contribute. 
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-1.94±0.23 


-I.59±0.2 


-1.06±0.17 


6 


-1.61±0.16 


-1.7±0.09 


-1.31±0.09 
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-1.53±0.37 


-1.39±0.2 


-1.15±0.16 
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-1.8±0.05 


-1.71±0.04 
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-1.63±0.63 
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-1.84±0.07 


-1.86±0.09 
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-1.75±0.12 
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-1.44±0.09 


49 


-1.99±0.38 


-1.65±0.16 


-1.36±0.12 
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-34.65±*** 


-23.7±0.91 


-23.47±0.32 
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-20.96±0.43 


-23.24±0.5 


-22.98±0.23 
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-21.79±0.83 


-22.55±0.69 


11 
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-22.19±0.37 
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-22.95±0.96 
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41 
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-23.38±1.5 
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49 


-31.66±*** 
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-23.72±0.46 
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57.51±*** 


24.66±17.2 


37.27±11.98 


6 


21.64±10.25 


41.53±20.52 


56.29±14.72 


8 


3.98±3.63 


3.07±1.92 


7.09±3.74 


11 


1±0.82 


1.11 ±0.97 


1.51±0.85 


13 


1.86±1.53 


12.93±4.75 


23.17±12.84 


18 


1.27±0.99 


107.34±*** 


14.61±6.25 


21 


0.94±*** 
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1.82±1.11 


22 


3.58±2.45 


4.22±1.64 


7.27±3.13 


25 


10.95±*** 


10.81±6.65 


14.78±6.95 


27 


7.19±5.75 
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14.81±8.72 


40 


5.2±2.95 


10.1±6.79 


12.87±5.47 


41 


12.99±*** 


12 21±*** 


5.06±4.05 


44 


95.01±*** 


9.45±4.25 


13.28±6.7 


49 




13.86±9.14 
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Table 2. Results of the Schechter function fitting of the LFs of the 14 CI 
galaxy clusters. For some clusters, the AI* values were not constrained by 
the fitting program and the errors of these unconstrained M* are starred. 
Also, the errors of the corresponding Lqpt values are starred, since the 
computation of Lqpt depends on both a and M*. 



To the right (fainter side) of this line, one or more of the clusters in 
the redshift bin drop out of the stacked data. 

To check whether the trend of a with redshift is robust against 
these two effects, we carried out tests on the stacked data, by pro- 
gressively removing the faintest magnitude bin in the stacked LFs 
and re-fitting. In general, we found no significant change in the fit- 
ted values of a (which changed only within their errors), or in the 
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0.92±0.11 


0.93±0.07 


0.90±0.08 


Lqpt.T 


0.95±0.06 


0.96±0.04 


0.97±0.03 


LopTjT.z (Partial) 


0.87±0.16 


0.89±0.11 


0.92±0.06 




0.82±0.21 


0.83±0.16 


0.86±0.10 


Lopy,z,T (Partial) 


0.28±0.44 


0.36±0.36 


0.61±0.22 




-0.08±0.30 


0.23±0.30 


0.64±0.21 


a,T 


0.01 ±0.29 


0.05±0.30 


0.54±0.25 


Q:,T,z (Partial) 


0.00±0.29 


-0.37±0.28 


-0.17±0.30 


a,z 


0.01±0.29 


0.44±0.27 


0.54±0.25 


a,z,T (Partial) 


0.20±0.30 


0.67±0.20 


0.65±0.21 


M*,T 


-0.32±0.43 


-0.57±0.31 


-0.48±0.26 


M*,T,z (Partial) 


-0.47±0.41 


-0.15±0.37 


-0.02±0.29 


M*,z 


-0.06±0.43 


-0.57±0.31 


-0.48±0.26 


M*,z,T (Partial) 


0.25±0.44 


0.00±0.36 


0.00±0.29 


Redshift-stacked 


a,z 


0.97±0.10 


0.88±0.29 


0.89±0.29 


a,z,^ (Partial) 


0.91±0.24 


0.65±0.58 


0.51 ±0.65 


M*,z 


-0.86±0.33 


-0.87±0.31 


-0.71±0.53 


M*,z,T (Partial) 


-0.21±0.67 


-0.58±0.62 


-0.33±0.68 


Temperature-stacked 


a,T 


0.10±0.64 


0.31±0.68 


0.75±0.49 


a,T,z (Partial) 


-0.44±0.67 


-0.85±0.35 


-0.45±0.67 


M*,T 


-0.46±0.67 


-0.26±0.67 


0.67±0.56 


M*,T,z (Partial) 


0.20±0.67 


-0.46±0.67 


0.25±0.67 



Table 3. Pearson's correlation coefficients (PC. Coeff.) of individual CI 
clusters, redshift-stacked clusters and temperature-stacked clusters for the 
three-filter set (g' , r',z'). 'X,Y,Z (Partial)' denotes partial correlation co- 
efficient of quantities X and Y with effects of quantity Z removed, to be 
compared with the line directly above it, where correlation coefficient of 
the same quantities X and Y is presented without partial analysis. 



a-z correlation when the LFs were truncated at the vertical dashed 
line, or when the LFs for all redshift bins were fitted to the same 
limiting absolute magnitude (which is set by the most distant sys- 
tems). There was one exception to this. The three clusters in the 
highest redshift band (clusters 1, 6 and 49) all have a values (al- 
beit with large errors) steeper that the shallow slope of -1.31 which 
fits to the stacked data in the g' band for this high redshift bin. As 
the faintest bins, to the right of the dashed line in the bottom plot 
of Figure [Tol are progressively removed, the fitted slope steepens. 
Hence the flat slope of -1.31 must be regarded as unsafe, and the 
very high a-z correlation in g' band, given in Table[3] is probably 
overestimated. Rather, we have a situation in all three photometric 
bands, where the faint-end slope is steeper al z < 0.2 than it is at 
higher redshift. 

To visualise the behaviour of the faint-end slope in terms of 
both redshift and colour, we plot the faint end of the fitted lumi- 
nosity functions for the three bands in Figure [3] using green for g' 
band, red for r' band and black for z' band. For this plot, we have 
divided the sample into three redshift bins: low (0.05-0.14), inter- 
mediate (0.26-0.32) and high (0.43-0.61), denoted by different line 
styles. All LFs have been renormalised to have = 1 at M=-19.5. 

The Figure shows how the faint end slope steepens towards 
low z. It also illustrates colour trends in a. At low redshift (solid 
lines), the slopes are very similar (though the curves are separated 
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Figure 3. The faint end of the fitted LFs of the C 1 sample grouped into 
three redshift bins: low (0.05-0.14, solid lines), intermediate (0.26-0.32. 
dashed lines) and high (0.43-0.61, dotted lines). Colours represent the filter 
bands: green for g', red for r' and black for z' . All LFs were normalised to 
have (j> = 1 aS. M=-19.5 for easy comparison. The faint end slopes become 
shallower with increasing redshifts. Also, at intermediate redshift (dashed 
lines), the slope shows a trend with colour, becoming steeper towards the 
blue. This colour trend largely vanishes at low redshifts (solid lines). 
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Figure 4. The faint-end slope of the fitted Schechter function in g', r' and 
z' bands for local clusters with redshift 0.05 to 0.14 and for intermediate 
redshift (0.29 < z < 0.32) clusters. 



due to their different values of M*), whilst at intermediate redshift 
(dashed lines), the slope shows a strong trend with colour. 

The values of a in Table|4]also show a trend with colour. The 
faint-end slope of ^-stacked clusters becomes steeper as we move 
from z' (red side) band to g' (blue side). This trend is very obvious 
in the second, third and fourth redshift bins (0.29 < z < 0.32) and 
much less obvious and maybe absent (within the errors) in the first 
hm(z < 0.14), see Figure|4]in which we plotted the values of a for 
the three bands for the lower- and intermediate redshift bins. The 
increase in the faint-end slope of the Schechter function in the bluer 
bands means that at the faint side of the colour-magnitude diagram 
the blue galaxies outnumber red ones. 

To explore this we produced K-corrected colour-magnitude di- 
agram (Figure Is} of g'2 — z'2 versus absolute rkron magnitude for 
0.29 < z < 0.32 (six clusters: 8,13,18,22,27 and 40) in which 
this trend is most obvious, and the same plot for the first redshift 



Figure 5. Colour-magnitude diagram: g' — z' versus r' (K-corrected) for 
low-redshift (z=0.05-0. 14) clusters and intermediate-redshift (z=0.29-0.32) 
clusters. 



bin, 0.05 < z < 0.14 in which no such trend is apparent. Figure|5] 
clearly demonstrate how the distribution of cluster galaxy colours 
changes from z ~ to z ~ 0.3. In the g' band the evolution of a is 
much stronger, especially after removing the effects of the temper- 
ature (partial correlation). These trends in a show that the fraction 
of blue faint galaxies at 2 ~ 0.3 was larger than it is now, and 
suggests that these galaxies have reddened and moved 

The Schechter function characteristic magnitude M* in the 
redshift-stacked clusters showed a negative correlation with red- 
shift. The correlation coefficients between M* and redshift are high 
but less significant than those between a and redshift. However, 
when the partial calculations were carried out, these coefficients 
dropped and became consistent with zero. Hence the trend in M* 
with z appears to be due to a selection effect: hotter clusters are 
more luminous, and so are more easily detected at high r edshift, and 
these brighter clusters also tend to have brighter M* ( IZandivare3 

(iooi)). 



3.4 Luminosity functions of 2 ~ 0.3 clusters 

Eight amongst the 14 CI clusters, more than half of our sample, 
lie within the narrow redshift range 0.26 to 0.32. These clusters 
are representative of low-mass clusters at intermediate redshifts - a 
population which dominates the XMM-LSS cluster dataset. Stack- 
ing these clusters together provides the best available composite LF 
for X-ray selected poor clusters at z ~ 0.3, which should be valu- 
able for future comparative studies. The LFs and their associated 
error contours are shown in Figures|6]|7]and[8] 

These clusters range in temperature from 1.3 to 2.8 keV. 
Schechter fits give a values -1.66±0.11, -1.50±0.05 and -1.36 ± 
0.05, and M* values -21.07±0.38, -22.21±0.22 and -22.83 ± 
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0.17, for the g', r' and z' bands respectively. Their faint-end slopes 
are shallower than the local clusters (z < 0.14) but steeper than 
higher-redshift (z > 0.43) ones. The colour trend of a is very ob- 
vious and seems to be a characteristic of 2 ~ 0.3 clusters com- 
pared to other cluster s in other redshift bins, as mentioned above. 
IWilman et alj i2005al) studied a sample of poor clusters at redshift 
z ~ 0.4 selected from the CN0C2 galaxy redshift survey. Compar- 
ing this optically selected sample with nearby clusters, they found 
that the fraction of passive galaxies declines with redshift, which is 
consistent with our finding of larger population of faint blue galax- 
ies at 2: ~ 0.3. However, these authors did not study the LF of their 

intermediate redshift g roups. 

In a recent study, iHarsono & ProDrij I2OO 9) presented com- 
posite LFs of six rich (T ~ 7-9 keV) clusters with redshifts rang- 
ing from 0.14 to 0.40 (averaging to 0.246) in the B,g,V,r,i and z 
bands. The LFs were well fitted by a single Schechter function 
with a values for g, r and z bands as follows: —1.31 ± 0.04, 
— 1.33 ± 0.03 and —1.45 ± 0.02 and the corresponding M* values 
were -20.94±0.17, -21.95±0.29 and -22.26±0.30. Their M* 
values are in reasonable agreement with ours, but their slopes are 
shallower, and show no trend with colour However, their data were 
limited to 20-40% of the area within r2oo, and they suggest that the 
lack of any upturn in the slope at faint magnitudes may be related 
to this - the extra faint galaxies responsible for the upturn being 
associated with a population infalling into clusters. In contrast, our 
data extend to l.Srsoo, which is approximately equal to r2oo. 



3.5 Temperature-stacked clusters 

The 14 CI clusters span a temperature range of 0.64 to 4.80 keV 
This was divided into five subranges: 0.64-1.00, 1.30-1.34, 1.60- 
2.20, 2.80-3.20 and 4.80-4.80, using the same criteria, discussed 
above, which were used for stacking into redshift bins. The highest 
temperature bin consists of only one cluster because after removing 
the three high-redshift clusters (2,5,29), the temperature difference 
between the two highest temperature clusters was too large to stack 
them together, and the LF of the highest temperature cluster, clus- 
ter 6 (T=4.80 keV) was of sufficient quality that it provides useful 
constraints on its own. The second highest temperature bin consists 
of two clusters and the rest have at least three clusters. 

The correlation coefficients of a with temperature are not high 
enough to establish any trend, especially when we take into consid- 
eration the reverse in sign of the coefficients after the partial corre- 
lation calculation, see Table[3] But in Table|4]the highest tempera- 
ture bin contains only one cluster (cluster 6) and the other bins show 
some indication that a increases (slope decreases) with tempera- 
ture, especially in bands r' and z' . Further investigation is needed 
to arrive at more conclusive results about the a tend with temper- 
ature. As to A'/*, the stacked results do not show any trend with 
temperature. The LFs of the temperature stacked data are shown in 
Figures[T3l[i4]and[T5] 

Some previous studies (see for example, iMiles et alj J2004h ) 
found that galaxy clusters with low X-ray luminosity (comparable 
to the coolest clusters in our CI sample) exhibit dips in their LFs. 
In our data, some of the temperature stacked LF plots ( 1131 [74l and 
\\5\ . especially those with high temperature (> 2.8 keV) showed 
signs of dips in the faint end of the LF. It can be hard to distinguish 
between scatter of the data points and a genuine dip in the LF. 

To test the genuineness of these dips we fitted a Schechter 
function minus a Gaussian function defined by three parameters 
(central magnitude, width and depth) to these temperature stacked 
LFs. The two fits, with and without the Gaussian were statistically 
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Figure 6. LFs of the 8 stacked CI clusters with redshift 0.2 to 0.4 and their 
associated la, la and Scr contours of confidence levels for a and M* in 
the g' band. 



compared using their x values, and an F-test applied to assess the 
significance of the improvement resulting from inclusion of the dip. 
In some cases the dip improved the fit at a confidence level of more 
than 90%. See, for example Figure|9] 

However, careful examination of the stacked LF and the indi- 
vidual clusters LFs in these cases suggested that the dip is produced 
by the stacking of clusters with different faintest magnitude limits, 
rather than lying within the magnitude range shared by all the com- 
bined clusters. This was found to be true for all stacked LFs that 
showed a statistical improvement in fit on inclusion of a Gaussian 
dip. We therefore conclude that our data show no evidence for real 
dips in the optical LFs of the CI clusters. 



4 DISCUSSION 

4.1 Faint-end slope of the luminosity function 

In this work we have studied the LFs of the individual clusters in 
the CI sample from XMM-LSS. A Schechter function provided a 
reasonable fit across most of the LF for most clusters, especially 
in the z band. But the bright end was poorly-fitted for nearly half 
of the sample (6 out of 14) and M* values were often not well 
constrained. The faint-end slope ranges are —1.03 < a < —2.1, 
-1.19 < a < -1.89 and -1.06 < a < -1.77 with averages 
-1.70 ±0.10, -1.64 ±0.04 and -1.43 ±0.03 for the g', r' and z' 
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Table 4. Results of the Schechter function fitting of the redshift-stacked and temperature-stacked clusters for the three-filter set (</', t' ,z'). 



passbands respectively. The mean faint-end slope, averaging over 
all the three filters, is aavg ~ —1.59 ± 0.05. 

Comparison of fitted Schechter parameters from different 
studies should take into account the passband, cluster redshift, and 
the procedure used in constructing the LF, including the methods 
used to determine the cluster membership and the background sub- 
traction, since all of these factors may affect the results and there- 
fore the accuracy of comparison. 

Previous studies of gala xy cluster LF s have found a wide 
range for a, from a ~ —1 dPaolillo et al.l (|2001)) to a ~ —2 
jPopesso et al.l ( I2006 )). but generally, LFs of clusters (both high- 
mass and low-mass systems) are found to have steeper sl opes than 
field g alaxy L Fs, which usua l ly spa n values a ~ —0.7 jLin et al] 
( Il996h ) to -1 jLovedav et al] jl995l) ). The values we obtain for a 
fall into the cluster LF range. The mass (M500) r ange of the CI 
clusters is 0.6-19 xIO^^Mq jPacaud et atl ilOOlh ) and this puts 
the CI sample in the lower-mass class of galaxy clusters (poor clus- 
ters and groups). This indicates that low-mass systems have almost 
the same range of faint-end LF slopes as more massive systems. 

Moreover, iGonzalez et alj j2006h studied LFs of galaxy clus- 
ters with a virial mass range 0.01 — 20 x 10^^ Mq and redshift 
0.03 < z < 0.06 and found slopes of -1.9 < a < -1.6 at 
the faint end (Mr > —18). This is consistent with our result for 
clusters with comparable redshift; clusters 11 and 21, which have 
estimated masses of 0.6 and 0.9 x lO" Mq and redshifts 0.05 and 
0.08, show faint-end slopes of -1.80 ± 0.05 and -1.89 ± 0.06 in 
the r' band, with magnitude limits of -14.5 and -15 respectively. 

The study of IPopesso et al.l ( l2005al) on X-ray selected rich 
clusters with z < 0.25 also gave a steep faint-end (Mg > —16) 
slope: -2.1 < a < -1.6 in the SDSS g band. CI clusters with 
redshifts < 0.26, namely clusters 11,21,25,41 and 44, have a g' 
ban d slope range of —2. 1 < a < —1.67, which agrees well 
with lPopesso et al] ilOOS^ . The redshift-stacked clusters with red- 
shift z < 0.32 (the first four redshift bins) gave a slope range 
of —1.7 9 < Q < —1.59 in the g' band which is also consis- 
tent with lPopesso et al] j2005ah . The C I clusters are low-mass sys- 
tems, whilst the iPopesso et 'Zl ( l2005iJ) systems are rich clusters. 
The steep faint end slopes seen in both indicate a larger fraction 



of dwarf galaxies in both groups and clusters, compared with the 
shallower LF slopes usually found in the field. 

However, as is the case with richer clusters, the results from 
different studies of the luminosity function of group galaxies arrive 
at different results. For instance. iMiles et al . (2004) derived a very 
flat (a ~ —1) Schechter slope for X-ray bright grou ps - though 
they f ound a faint upturn in X-ray dim systems - and lZandivarezI 
(I2OO6I) derived sim ilarly low faint end slopes for SDSS groups. 
iMiles et alj ( |2004|) used photometric data of X-ray selected sys- 
tems and used all galaxies with B-R < 1.7 from the regions outside 
a radius of i?5oo from the centre of the group as the background 
for subtraction, whilst Izandivarezi j2006h used spectroscopic data 
for mem bership determination f or their friends-of-friends selected 
clusters. iRobotham et al.l ( I2OO6I) extracted LFs for 2PIGG groups, 
derived from the 2dF galaxy redshift survey, and obtained good 
fits with Schechter functions, with faint end slopes which increased 
from a ~ — 1 for red galaxies to a ~ — 1.5 for blue galaxies. These 
discrepancies in the faint-end slope from different studies could 
arise from a variety of causes: different cluster selection methods 
(X-ray selected clusters in our case), spectroscopic or photometric 
selection of cluster galaxies, different galaxy background subtrac- 
tion techniques (see discussion in section 4.3), and possibly be- 
cause different clusters have different faint-end slopes depending 
on their large-scale environment, which will affect the incidence of 
inf ailing galaxies. 



4.2 The absence of the upturn in the faint end of LFs 

Both lPopesso et af] ( l2005al) and I Gonzalez et ai] j2006l) reported an 
upturn at the faint end of their stacked LF, and required a sum of 
two Schechter f unctions, rather than a single Schechter, to obtain 
reasonable fits. iPopesso et al.l j2005ah locat ed the upturn at -16 in 
the g' band, and -18.5 in z'; the upturn of iGonzalez et alj ^20061) 
was found at a similar magnitude: -18 in the r' band. In our sam- 
ple, only the LFs for clusters 11,21 and 41 e xtend to the faint mag- 
nitudes in which IPopesso et ei\ ( l2005al) and IGonzalez et af] j2006h 
found their upturns. The composite LF for these systems is the first 
in the redshift-stacked LFs, see F iguresllOlfTTIandfni Althou gh our 

Syndic " 



results agree with lPopesso et al] ( l2005af) and Gonzalez et alj 1 20061) 
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Figure 7. LFs of the 8 stacked CI clusters with redshift 0.2 to 0.4 and their 
associated la, 2a and 3a contours of confidence levels for a and M* in 
the r' band. 



Figure 8. LFs of the 8 stacked CI clusters with redshift 0.2 to 0.4 and their 
associated la, la and Scr contours of confidence levels for a and M* in 
the z' band. 



regarding the steep values of the faint-end slope, we do not find any 
evidence for a departure from a simple power law at the faint end. 

Other studies gave steep slopes at the faint end of cluster 
LFs but wi t hout evidence of sudden upturns, see for example 
iDurret et alj j2002h . lOarilli et all fl999l) studied composite LFs of 
65 clusters ranging in redshift from 0.05 to .25 and did not find 
upturns in their com posite LFs. lPopesso et all (2005a) argued that 
iGariUi et alj h999l) did not see this upturn in their stacked LF be- 
cause they used a weighting for the individual LFs which depends 
strongly on the cluster magnitude limit, such that clusters with 
fainter magnitude limits, which contribute to the faint magnitude 
bins of the stacked LF, were heavily down-weighted. In our sam- 
ple, only the LFs for clus ters 11,21 and 41 ext end to the f aint mag- 
nitude region in which iPopesso et al. I j2005ah and Gonzalez et al.l 
( l2006h found their upturn. The composite LF for these systems is 
the first in the redshift- stacked LFs, see Figures fTOlll ll andll2l We 
did not apply any weighting method that depends on the magnitude 
limit and although faint-end slopes are steep in all three bands for 
the stacked LF of clus ters 11,21 and 41, the y lac k any upturn at 
the lo cations found by IPopesso et alj i2005ah and lOonzalez et al.l 
(12006^. Furthermore, individual LFs of these three clusters do not 
show any obvious upturn in the faint-end part of the LF that can 
be distinguished from the scatter of the data relative to the fitted 
Sch echter function . 

IPopesso etal.l ( l2006l) decomposed their LF by galaxy type and 
showed that the late-type galaxies LF was well fitted by a single 



Schechter function with a steep slope (a = 2.0 ), while the early- 
type galaxies LF could not be fitted by a single Schechter function, 
and a composite of two Schechter functions was needed, such that 
the faint-end upturn of the global cluster LF was due to the early- 
type cluster galaxi es. This sugges t s one way of reconciling our re- 
sults with those o f IPopesso et aP j2006h . If in our poorer clusters 
late-type galaxies outnumber early-types in the intermediate and 
faint magnitude ranges then the LF would be steep and without any 
upturns. This needs to be further investigated by studying the early- 
type and the late-type LF separately. Anoth er possibility for th e dif- 
ference between our results and those of iPopesso et alj ( 12006 ) lies 
in the techniques used to remove non-cluster galaxies, as we dis- 
cuss in the next section. 

4.3 Membership determination methods: Effects on a 

The steepness of the faint end of the luminosity function reflects 
the number of dwarf galaxies within a cluster. Estimates of this 
number are very sensitive to the method used to estimate and re- 
move the contribution of background and foreground galaxies be- 
fore constructing the cluster luminosity function. 

Rines & Geller (2008) compared methods of membership de- 
termination based on spectroscopic data and on photometric data 
(which we used) with regard to the resulting LF. They highlighted 
the advantages of spectroscopic identification of cluster members. 
Where automated photometric methods are used, they found, for 
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Figure 9. LF of the stacked clusters 1 and 27 (fourth temperature bin: 
2.80 < T < 3.20) with plot of the Icr, 2(t and 3cr confidence contours 
in the a-M* plane. The line is the fit of a Schechter function plus a Gaus- 
sian dip. The fitted dip position is — 19.9±0.1 (z' filter). The vertical dotted 
line marks the faintest magnitude at which both stacked clusters contribute. 



4.4 Origin of the faint galaxies 

Our results indicate that larger numbers of faint galaxies exist in 
cluster environments than in the field. It is not straightforward to 
understand this result, since various dynamical processes that can 
destroy dwarfs act more effectively in dense environments. Several 
idea s have been proposed to explain the excess of dwarfs in clus- 
ters. i Babul & Reel j 19921) argued that a primordial population of 
dwarf galaxies is preserved in high-pressure environments, whilst 
it fades away in low-pressure regions. Alternatively, dwarfs could 
be formed by galaxies that fell into clusters from the surrounding 
field and were morphologically transformed. The transformation 
mechanism could be tidal fragmentation or so-called harassment 
of inf ailing late-t ype spiral ga l axies by the cluster potential or by 
close encounters (Moore et al. ( 1996) ) or ram pres sure stripping of 
dwa rf irregular ga l axies (e.g., van Zee et al. 

iBoselli et al. showed that both simulations and obser- 

vations are consistent regarding the scenario of recent accretion and 
transformation of low-luminosity star-forming galaxies in the Virgo 
cluster into quiescent dwarfs due to ram pressure gas stripping and 
galaxy starvation. They also showed that this process of transforma- 
tion results in galaxies with structural and spectrophotometric prop- 
erties similar to those of dwarf ellipticals. If the whole star-forming 
dwarf galaxy population dominating the faint end of the field lu- 
minosity function were accreted, it could be totally transformed by 
the cluster environment into dwarf ellipticals on timescales as short 
as 2 Gyr. These vigorous forces acting in cluster environments may 
explain the steepness of LPs faint-end slopes in nearby clusters. 



example, that many large galaxies, especially those with low sur- 
face brightnesses, may be detected as many small separate objects, 
and warned that if these pieces of galaxies are not removed, they 
can produce an artificial excess of faint galaxies in cluster fields. 

However, we have to emphasise that although spectroscopic 
data can give precise information on the cluster membership, their 
use to study cluster LPs is limited to relatively nearby clusters, 
since for higher-redshift clust ers, spectroscopy is feasible only for 
the bright cluster galaxies. iBoue et al. I (2008) used deep multi- 
colour photometry to study the LP of A496, using colour selec- 
tion to reduce contamination by red background galaxies, and did 
not find the large fraction of dwarf gal axies (a = 2 .0) inferred 
by some other authors, including Popessoetah I l l2006l) . They sug- 
gested that this excess of dwarf galaxies in some studies might arise 
from inadequate removal of background, due to use of inadequate 
(or no) colour cuts. They claimed that the red sequence used by 
IPopesso et al.l ( |2006|) was polluted by field galaxies because they 
used u* — r' vs i' in their CMD which Boue et al. (2008) showed 
was not efficient in rejecting background galaxies. In our study, we 
did not use u* — r' vs i' to define the colour cuts. Instead, we used 
u* — g' vs g' for the g' band, g' — r' vs r' for the r' band and 
i' — z' vs z' for the z' band. Moreover, our method of field LP 
subtraction is based on global background LP constructed by using 
the whole 1° x 1° field of the cluster. Therefore, we don't see any 
obvious reason why we might have contaminated the red sequence 
with field galaxies in such away as to give a false steep faint-end 
slope. 



4.5 The evolution of a 

Our results show an evolutio nary trend of the faint-end slope, a, in 
all bands used: g', r' and z'. iLiu et al.l ( 1200^ examined the faint- 
end slope of the V-band LP of field galaxies with redshifts 2 < 0.5 
and found that it becomes shallower with increasing redshift: their 
a changed from — L24 for the lowest redshift bin 0.02 < z < 0.1 
to —1.12 for the highest redshift bin 0.4 < z < 0.5. In clusters, 
a recent study by ILu et al.l ( |2009|) of an optically selected cluster 
sample found steepening of the faint end with decreasing redshift 
since z ~ 0.2, and that the relative number of red-sequence dwarf 
galaxies had increased by a factor of ^ 3. 

It is pos sible that this LP slope tren d with redshift is linked to 
the finding o arsono feProprislHoOTh that the 'upturn' in the LP 
faint end (i.e. the excess of galaxies above a single Schechter func- 
tion) is found only in low redshift clusters. They attributed this to 
the recent infall of star-forming field galaxies or the whittling down 
of formerly more massive objects. The impact of rec ent infall of 

! galaxies into clusters is also supported by the work of iLisker et al.l 
20071) ■ who showed that dPs in the Virgo cluster fall into two ma- 
jor morphological subclasses: a) dPs with blue centres, thick disks 
or features reminsiscent of late-type galaxies, such as spiral arms or 
bars; this class showed no central clustering, suggesting that they 
are an unrelaxed population formed from infalling galaxies. The 
second subclass is b) nucleated dPs - a fairly relaxed population of 
spheroidal galaxies indicating that they have r esided in the cluster 
for a long time, or were formed along with it. iLisker et alj j2007h 
also pointed to other studies deriving similar results (see references 
therein), indicating that this subclassification is not specific to the 
Virgo cluster. 
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4.6 Colour trends 

The faint-end slopes, a, of the redshift- stacked groups are steeper 
in bluer bands in almost all redshift bins. However, this trend is 
significant (> la) only for the redshift range 0.29 < z < 0.32. 
The redshift bin in which this effect seems to be absent is the first 
bin: 0.05 < z < 0.14. This suggests that the fraction of faint 
blue galaxies in clusters of redshift z ~ 0.3 are higher than in lo- 
cal systems. Figure |5] further illustrates this and it also shows that 
these blue faint galaxies appear to have reddened and moved up- 
wards in th e colour-magn i tude dia gram. This is consistent with the 
findings of IWilman et~S] ( l2005ij) who compared the fractions of 
passive (red and quiescent) and blue star-forming galaxies in clus- 
ter at 0.3 < z < 0.55 with nearby (z ~ 0) clusters. They found 
that the fraction of passive galaxies declined strongly with redshift 
to at least z ~ 0.45. These results are also consistent with the well- 
known Butcher-Oemler effect in clusters and support the idea that 
dense environments are responsible for galaxy transformation from 
blue to red becaus e these trends are less obvious in field environ- 
ments, see Wi lman et alj ( l2005bl) . 

Our result is also consistent with lYee et al.l ( |2005|) who stud- 
ied the colours of galaxies as a function of luminosity and environ- 
ment using the Red Sequence Cluster Survey and the SDSS. They 
found a higher incidence of faint to moderate luminosity galaxies in 
high density environments at z > 0.2 compared to lower redshifts 
and lower density environments. They interpreted this as arising 
from the shut-down of star formation in low mass galaxies within 
clusters at z < 0. 3, in contrast to the situation in the field (c.f. 
[Balogh etal]( l2004h ). 

The fact that such transformations are observed in low-mass 
clusters like our CI sample, as well as in richer clusters, favours 
mechanisms for suppression of star formation which operate in 
shallower potential wells, such as strangulation, tidal interactions 
and galaxy mergers, rather than ram pressure stripping, which is ef- 
fective mostly in rich environments with high velocity dispersions. 

4.7 Correlation between global properties of clusters 

The optical luminosity is a good indicator of cluster richness, and 
hence should be clos ely related to cluster mass, velocity disper- 
sion and temperature jPopesso et al] j2005bh ). Assuming that clus- 
ter mass is directly proportional to the optical light (i.e. M/ Lqpt 
is constant), that the ICM is in hydrostatic equilibrium and that X- 
ray luminos i ty Lx scales with gas temperature T as Lx oc T'^ 
jXue & Wul bOOd) ). it is expected that Lqpt oc T^'^, and that 
LoPT oc Lx^. 

Our scaling results for Lqpt with Lx and T mostly agree 
well with these expectations. For the Lqpt — Lx relation, the 
logarithmic slopes are 0.47 ± 0.07 {g' band), 0.43 ± 0.08 (r' band) 
and 0.50 ± 0.07 (z' band). While for the Lqpt - T relation, we 
have 1.57 ± 0.17 (g' band), 1.51 ± 0.17 (r' band) and 1.79 ± 0.12 
(z' handy 

Pope sso et all j2004l) found 0.38 ± 0.02 for the Lqpt - Lx 
relation and 1.12 ± 0.08 for the Lqpt — T relation in the z 
SDSS band within a cluster radius of 0.5 Mpc (chosen to minimise 
the scatter in their scaling relations). The systems they used for 
their analysis, the RASS-SDSS sample, were X-ray selected, and 
ranged from low-mass systems of to massive clusters 

of 10^^ Mq, over a redshift range from 0.002 to 0.45. Their loga- 
rithmic slope value for the Lqpt — Lx relation is not inconsis- 
tent with our value (within the errors), however, their Lo pt ^ T 
value is lower than ours. They attributed the departure of their re- 



sults from the expected values to a breakdown in the assumption of 
constant mass-to-light ratio. More precisely, they argued that if the 
assumption of hydrostatic equilibrium was retained, their results 
wo uld be in a good agr eement with M/L oc i,o.'i3±o.o'i ^ found 
bv lOirardietaU ( |2002|) . However, we note that extracting Lqpt 
within a fixed metric radius, will include a smaller fraction of the 
virial radius for higher mass clusters. Hence it should be no surprise 
if the Lo PT ~ T relation is flattened below the expected slope of 
1.5 for self-similar clusters. U sing clusters from the RASS-SDSS 
sample, IPopesso et all bOOSbh calculated Lqpt within rsoo and 
r2oo- Their r2oo results were (we used 1.5 x rsoo): 0.57 ± 0.03 
(g band), 0.58 ± 0.03 (r band) and 0.58 ± 0.03 (z band) for the 
LoPT - Lx relation and 1.62 ± 0.10 {g' band), 1.64 ± 0.09 (r 
band) and 1.62 ± 0.10 (z' band) for Lqp t — T relation . Thes e 
are in better agreement with our results than lPopesso et al.] ( l2004h . 
and demonstrate the importance of the radius used to estimate the 

LoPT- 



5 CONCLUSIONS 

We have studied the luminosity functions of 14 Class 1 (CI) XMM- 
LSS galaxy clusters in three CFHTLS MegaCam bands: g', r' and 
z' . The X-ray selected clusters have masses ranging from 0.6 to 19 
X 10^^ M0, a redshift range of 0.05 to 0.61, and ICM temperature 
range of 0.64 to 4.80 keV. We used colour-magnitude lower and 
upper cuts to reduce contamination by cluster non-members, and 
performed background subtraction using the 1° x 1° field of view 
in which the cluster lies. K-corrected luminosity functions of galax- 
ies within 1.5 X r2oo were constructed for each cluster and fitted 
with a Schechter function. Total optical luminosities of the indi- 
vidual clusters were also computed by integrating over the fitted 
Schechter functions. The individual LFs were also stacked together 
into 5 different redshift and temperature bins. The main findings 
are: 

• A Schechter function provides a good fit across most of the 
LF for the majority of clusters in our sample. The value of a range 
from —1.03 to —2.1, but no evidence is found for upturns at the 
faint end of the Schechter function, even in the lowest redshift sys- 
tems, for which our LFs extend well into the dwarf regime. 

• M* ranges from —19.66 to —24.31. However, for many 
(nearly a third) of the clusters' M* values are not well-constrained. 

• The redshift-stacked LFs confirm that a becomes shallower 
with increasing redshift. The value of a is —1.75 ± 0.02 at low 
redshift (0.05-0.14), flattening to -1.22 ± 0.06 at high redshift 
(0.43-0.61) in the z' band. Similar trends are present in the other 
two bands. 

• a, also steepens significantly from the red (2') to the blue {g') 
band for clusters at redshift ~ 0.3. This effect is not present in 
our local clusters (z ~ 0), suggesting reddening of the faint blue 
galaxies from z ~ 0.3 to 2: ~ 0. 

• After removing the effects of redshift (correcting for the 
Malmquist effect), the temperature-stacked LFs do not exhibit any 
strong evidence for trends of the Schechter parameters with ICM 
temperature. 

• Total optical luminosities for our sample range from 1.0 to 
56.3 xIO^^Lq, and correlate strongly with X-ray luminosity. The 
logarithmic slopes of the Lqpt — Lx relation are 0.47 ± 0.07, 
0.43±0.08 and 0.50±0.07 for the g' , r' and z' bands respectively. 

• Also, LoPT correlate strongly with the X-ray gas temper- 
ature, T. The logarithmic slopes of the Lqpt — T relation are 
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1.57 ± 0.17, 1.51 ± 0.17 and 1.79 ± 0.12 for the g\ r' and z' 
bands respectively. 

• The slopes of the Lo pt — Lx and Lo pt — T relations are 
consistent with the established, non-self-similar, cluster Lx —T re- 
lation and constant mass-to-light ratio, except for the z' band value 
of the LopT — T relation which is higher than the expected value 
(1.5) by ~ 0.3. 

• Some of our stacked LFs show dips, but these appear to be 
artefacts arising where clusters with different faintest magnitude 
limits are stacked together. We therefore we conclude there is no 
evidence for real dips in the optical LFs of the CI clusters. 
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Figure 10. LFs of the stacked clusters for 5 redshift ranges and the associated la, 2cr and 3cr contours for g' band. All stacked clusters contributed to all 
magnitude bins are at the left side (brighter side) of the vertical dotted line which is at the faintest common magnitude bin of the clusters. Whereas at the right 
side (fainter side) of it, some clusters did not have data in some magnitude bins because they already reached their faintest magnitude limit. 
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Figure 11. LFs of the stacked clusters for 5 redshift ranges and the associated Ict, 2a and 3(t contours for the r' band. The vertical dotted line is at the faintest 
common magnitude value of all stacked clusters. 
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Figure 12. LFs of the stacked clusters for 5 redshift ranges and the associated Icr, 2cr and 3o" contours for the z' band. The vertical dotted line is at the faintest 
common magnitude value of all stacked clusters. 
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Figure 13. LFs of the stacked clusters for 5 temperature ranges and the associated Icr, 2cr and 3(T contours for the g' band. The vertical dotted line is at the 
faintest common magnitude value of all stacked clusters. 
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Figure 14. LFs of the stacked clusters for 5 temperature ranges and the associated Icr, 2cr and 3cr contours for the r' band. The vertical dotted line is at the 
faintest common magnitude value of all stacked clusters. 
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Figure 15. LFs of the stacked clusters for 5 temperature ranges and the associated Icr, 2rj and 3(t contours for the z' zand. The vertical dotted line is at the 
faintest common magnitude value of all stacked clusters. 
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Figure Al. LFs of the 14 individual CI clusters and contours of the well-fitted clusters for the r' band. Contours plots of the Itr, 2cr and 3cr confidence levels 
of a and M* are placed next to their associated LF. Clusters with failed constrained M* (and no contours) were placed at the bottom. 
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